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ABSTRACT 

Lamb-wave  testing  for  structural  health  monitoring  is  complicated  by  the  dispersion  nature  of  the  wave 
modes.  The  dispersion  effect  will  result  in  a  propagated  wave  with  longer  time  duration,  deformed 
envelop  shape  as  compared  to  its  excitation  counterpart,  and  hard  to  be  interpreted.  This  paper  first 
reviews  the  dispersion  compensation  and  removal  algorithms.  Second,  it  compares  these  two  methods 
by  applying  them  to  two  widely  used  low-frequency  Lamb  wave  modes:  SO  and  AO.  Numerical 
simulations  are  compared  in  parallel  with  experimental  results.  Finally,  the  dispersion  compensation 
algorithm  is  applied  to  1-D  PWAS  phased  array  and  demonstrated  to  improve  the  phase  array’s  spatial 
resolution. 
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1  INTRODUCTION 

Lamb- wave  testing  for  SHM  is  complicated  by  the  dispersion  nature  of  the  wave  modes. 
Dispersive  waves  have  frequency-dependant  propagation  characteristics  (Sachse  and  Pao,  1978).  Hence, 
even  narrowband  Lamb  wave  excitations,  such  as  tone  bursts,  will  disperse  as  they  propagate  through 
structures.  The  dispersion  effect  will  result  in  a  propagated  wave  with  longer  time  duration  and  deformed 
envelop  shape  as  compared  to  its  excitation  counterpart.  This  deteriorates  the  wave  spatial  resolution  and 
makes  it  hard  to  interpret  the  experimental  data. 

This  dispersion  issue  can  be  addressed,  for  example,  by  using  Lamb  wave  tuning  technique 
(Giurgiutiu  2005,  Santoni  et  al.  2007),  which  utilizes  narrowband  excitation  and  tune  excitation  frequency 
until  a  quasi  non-dispersive  Lamb  wave  mode  is  obtained.  Alternatively,  it  can  be  addressed  by  using 
Lamb  wave  dispersion  compensation  algorithm  proposed  by  Wilcox  (2001,  2003).  This  algorithm  makes 
use  of  a  priori  knowledge  of  the  dispersion  characteristics  of  a  guided  wave  mode  and  performs  signal 
processing  algorithm  to  map  signals  from  the  time  domain  to  the  spatial  domain  and  reverse  the 
dispersion  process.  The  basic  idea  underlying  dispersion  compensation  is  very  similar  to  the  time  reversal 
procedure  presented  by  Xu  and  Giurgiutiu  (2007),  in  which  a  signal  recorded  at  the  receiver  side  is  time 
reversed  and  propagated  back  to  the  source,  the  signal  is  then  compressed  to  its  original  shape  Instead  of 
physically  propagating  back  the  received  signal  to  its  source,  dispersion  compensation  method 
analytically  maps  the  received  signal  back  to  its  source  to  retrieve  its  original  shape  and  location  through 
signal  processing  approach.  More  recently,  Liu  and  Yuan  (2009)  proposed  a  dispersion  removal 
procedure  based  on  the  wavenumber  linear  Taylor  expansion.  By  transforming  a  dispersed  signal  to 
frequency  domain  and  interpolating  the  signal  at  wavenumber  values  satisfying  the  linear  relation,  the 
original  shape  of  the  signal  can  be  recovered. 

This  paper  first  reviews  the  dispersion  compensation  and  removal  algorithms.  Second,  it  compares 
these  two  methods  by  applying  them  to  two  widely  used  low-frequency  Lamb  wave  modes:  SO  and  AO. 
Numerical  simulations  are  compared  in  parallel  with  experimental  results.  Finally,  the  dispersion 
compensation  algorithm  is  applied  to  1-D  PWAS  phased  array  and  demonstrated  to  improve  the  phase 
array’s  spatial  resolution. 
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2  THEORY  OF  DISPERSION  COMPENSATION  AND  REMOVAL 
2.1  Dispersed  wave  simulation 

Figure  1  illustrates  the  use  of  PWAS  in  pulse-echo  mode  to  detect  crack  damage  in  a  thin- wall 
structure.  PWAS  transducer  was  bonded  to  the  structure  to  achieve  direct  transduction  of  electric  energy 
into  elastic  energy  and  vice-versa.  In  Figure  1,  the  compact  wavey(^)  denotes  the  forward  Lamb  wave 
excitation;  the  elongated  wave  g(t)  denotes  the  dispersed  Lamb  wave  reflected  from  the  crack.  The 
propagatedy(0  waveform  at  arbitrary  spatial  and  temporal  location  (x,  t)  can  be  predicted  as 

u{x,t)  =  —  r  =  -L  f”  F(co)e^^‘^*-’^^do)  (1) 


Figure  1 


principles  of  pulse-echo  method  to  detect  a  crack  near  PWAS  transducers  on  a  thin-wall 
structure 


Assuming  constant  reflection  coefficient  A,  and  xq/2  distance  separation  between  the  PWAS  and  the  crack; 
the  reflected  waveform  after  xq  propagation  distance  can  be  predicted  by 

g{t)  =  u{x,t)L  =^\  (2) 

Where,  F{cd)  is  the  Fourier  transform  of  f{i),  k  is  the  wavenumber.  Eq.  (2)  enables  us  to  simulate  a 
dispersed  wave  after  a  certain  propagation  distance. 

2.2  Dispersion  compensation  algorithm 


If  we  propagate  backward  the  g{t)  to  its  source  location  (i.e.,  set . 
direction  (i.e.,  set  x  =  -x),  we  have  dispersion  compensated  waveform  h{x)  as 


uui  111  icvci&cu  piupa^auuii 


h(x)  =  u(x,t) 


^=0  271^-^  t=o 


(3) 


where  G{co)  is  the  Fourier  transform  of  g{t).  Note  that  Eq.  (3)  is  the  fundamental  dispersion  compensation 
equation.  It  maps  the  signal  from  time  domain  to  spatial  domain  and  reverses  the  dispersion  process. 
Recall  the  definitions  of  group  velocity  {cd)  and  phase  velocity  {cd)  ,  we  have 


dcD-c^^{cci)dk 

(o  =  Cpi,ico)k 

Substitute  (4)  into  (3),  we  have 

h{x)  =  —  r  H{k)eJ^dk 

2n^-^ 


(4) 

(5) 


where,  H{K)  =  G{c6)c^  {co) ,  co  =  co{k) .  The  wavenumber  k  can  be  thought  of  as  a  spatial  frequency 

variable.  Thus,  x  and  k  are  dual  variables  for  transforming  between  spatial  and  wavenumber  domains,  in 
the  same  sense  as  t  and  co  are  dual  variables  for  transforming  between  time  and  frequency  domains. 
Inverse  Fourier  transform  methods  (e.g.,  IFFT)  can  be  used  to  calculate  h{x)  in  Eq.  (5).  However,  besides 
interpolating  G{co)  and  (co) ,  careful  setup  of  the  variables  in  spatial/wavenumber  domains  w.r.t.  those 

in  time/frequency  domains  is  needed  to  ensure  the  calculation  accuracy. 


2.3  Dispersion  removal  algorithm 


The  Lamb  wave  nonlinear  wavenumber  k{co)  causes  phase  distortions  and  results  in  a  dispersed 
waveform.  If  the  wavenumber  is  in  a  linear  relation  w.r.t.  co ,  there  will  be  no  dispersion  in  a  propagated 
waveform.  Based  on  this  concept,  the  dispersion  removal  algorithm  maps  the  signal  to  wavenumber 
domain  where  linear  wavenumber  relation  is  satisfied  to  removal  the  dispersion  effect.  The  linear 
wavenumber  relation  is  approximated  by  using  the  Taylor  expansion  of  the  wavenumber  k{(D)  at  the 
excitation  central  frequency  co^  up  to  the  first  order  as 

Koj)  ~  kin  {co)^k{cok\  +  k  (®o  )(a>  -  ^  (6) 

Cgr(®o) 

Therefore,  to  remove  the  dispersion  in  a  dispersed  wave,  such  as  g(t),  the  procedure  can  be  summarized  as 
follows: 

•  Apply  Fourier  transform  to  g(t)  and  get  G(k) ,  where  G(k)  =  G(co)  since  ^  is  a  function  of  co ; 

•  Use  Eq.  (6)  to  calculate  linear  wavenumber  values 

•  Interpolate  [A:,G(A:)]  pair  at  ki-^(co)  to  get  ; 

•  Apply  inverse  Fourier  transform  to  G[ki-^ (^y)]  to  get  a  wave  with  the  dispersion  removed. 

As  compared  to  the  dispersion  compensation  algorithm,  the  dispersion  removal  algorithm  needs  less 
computation  efforts. 

3  DISPERSION  COMPENSATION  VERSUS  DISPERSION  REMOVAL 

Dispersion  compensation  and  dispersion  removal  algorithms  are  compared  in  this  section.  The 
comparison  is  based  on  both  numerical  and  experimental  results  for  recovering  two  dispersed  SO  and  AO 
Lamb  waves 

3.1  Numerical  simulation 

Numerical  simulations  comparison  to  recover  the  two  widely  used  So  and  Aq  modes  on  a  3 -mm  and 
1-mm  aluminum  plates  (p  =  2780  kg/m^,  E  =  72.4x10^  Pa)  were  performed.  To  maximize  the  dispersion 
effect  of  one  mode  while  suppressing  the  other  mode,  tone  burst  excitations  were  selected  to  center  at  350 
MHz  for  So  mode  on  the  3 -mm  plate  with  group  velocity  Cgr  =  5380  m/s  and  at  36  kHz  for  Ao  mode  on  the 
1-mm  plate  with  group  velocity  =  1 163  m/s 

Simulation  results  of  So  mode  using  dispersion  compensation  and  removal  algorithms  are  plotted  in 
Figure  2.  Figure  2a  shows  a  Hanning  windowed  tone  burst  excitation  centered  at  350  kHz,  which  is  used 
to  excite  So  mode  Lamb  wave.  Figure  2b  illustrates  the  dispersed  So  mode  waveform  after  propagation 
distance  of  300  mm,  simulated  by  Eq.  (2).  Because  of  dispersion,  the  waveform  changes  its  envelope 
shape  as  it  propagates  through  the  structure.  Figure  2c  shows  the  waveform  recovered  by  dispersion 
compensation  algorithm  using  Eq.  (5).  Spatial  resolution  of  the  wave  packet  is  largely  improved  as 
compared  to  its  dispersed  version.  Meanwhile,  the  waveform  is  exactly  sitting  at  the  expected  spatial 
location  x  =  300  mm.  Figure  2d  shows  the  recovered  waveform  by  dispersion  compensation  algorithm 
and  Figure  2e  shows  the  recovered  waveform  by  dispersion  removal  algorithm  in  time  domain.  Both  of 
the  recovered  waveforms  are  very  close  to  their  original  tone  burst  excitation.  A  close  examination  of 
these  two  recovered  waveforms  shows  that  dispersion  removal  algorithm  show  slightly  better 
performance. 

Similar  simulation  results,  as  shown  in  Figure  3,  were  observed  when  applying  the  dispersion 
compensation  algorithm  to  a  dispersed  36  kHz  Ao  mode  wave  on  a  1-mm  aluminum  plate.  The  spatial 
resolution  of  the  dispersed  Aq  mode  wave  was  improved  after  mapping  it  from  time  domain  to  spatial 
domain  using  dispersion  compensation  algorithm.  However,  the  compensated  waveform  in  time  domain 
(Figure  3d)  seems  to  possess  higher  frequency  components  as  compared  to  its  original  tone  burst 
excitation.  This  may  due  to  the  artifacts  introduced  by  mapping  procedure  (e.g.,  interpolation  method)  in 


dispersion  compensation  algorithm.  In  contrast,  the  recovered  waveform  by  dispersion  removal  algorithm 
is  very  close  to  the  shape  of  its  original  tone  burst  excitation,  as  shown  in  Figure  3d. 
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Figure  2  Numerical  simulation  of  dispersion  compensation  of  350  kHz  SO  mode  on  a  3-mm  aluminum 
plate:  (a)  3.5-count  Hanning  windowed  tone  burst  center  at  350  kHz;  (b)  dispersed  SO  mode 
wave  after  x  =  300  mm  propagation  distance,  simulated  by  Eq.  (2);  (c)  recovered  SO  mode 
wave  in  spatial  domain  by  dispersion  compensation  algorithm;  (d)  recovered  SO  mode  wave 
in  time  domain  by  dispersion  compensation  algorithm;  (e)  recovered  SO  mode  wave  in  time 
domain  by  dispersion  removal  algorithm 
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Figure  3  Numerical  simulation  of  dispersion  compensation  of  36  kHz  AO  mode  on  a  1-mm  aluminum 
plate:  (a)  3-count  Hanning  windowed  tone  burst  center  at  36  kHz;  (b)  dispersed  AO  mode 
wave  after  x  =  400  mm  propagation  distance,  simulated  by  Eq.  (2);  (c)  recovered  AO  mode 
wave  in  spatial  domain  by  dispersion  compensation  algorithm;  (d)  recovered  AO  mode  wave 
in  time  domain  by  dispersion  compensation  algorithm;  (e)  recovered  AO  mode  wave  in  time 
domain  by  dispersion  removal  algorithm 


3.2  Experimental  verification 
3.2.1  Experimental  setup 

Figure  4  shows  Lamb  wave  dispersion  compensation  and  removal  experimental  setup  using 
permanently  bonded  PWAS  transducers.  It  consists  of  a  HP33120  function  generator,  a  Tektronix  5430B 
oscilloscope  and  a  PC  (Figure  4a).  Two  specimens  were  used:  one  is  a  1524mmxl524mmxlmm 
aluminum  plate  bonded  with  two  round  7-mm  diameter  PWAS,  400  mm  apart  (Figure  4b);  the  other  one 
is  a  1060mmx300mmx3mm  aluminum  plate  bond  with  two  7-mm  square  PWAS,  300  mm  apart  (Figure 
4c).  To  eliminate  boundary  reflection  interference,  modeling  clay  was  put  around  specimen  edges.  Instead 
of  using  one  transmitter  and  one  reflector  working  a  pulse-echo  mode,  a  pair  of  transducers  working  in 
pitch-catch  mode  was  used  for  better  SNR.  All  the  signals  from  receivers  were  recorded  by  oscilloscope 
and  saved  in  Excel  data  format. 


7mm  square  PWAS 


Figure  4  Dispersion  compensation  experimental  setup  and  specimens:  (a)  dispersion  compensation 
experimental  setup;  (b)  1524mmx1524mmxlmm  2024-T3  aluminum  plate  bonded  with  two 
round  7-mm  PWAS,  400  mm  apart;  (c)  1060mmx300mmx3mm  2024-T3  aluminum  plate 
bond  with  two  7-mm  square  PWAS,  300  mm  apart 

3.2.2  Experimental  results 

Before  applying  dispersion  compensation  and  removal  algorithms,  the  recorded  waveforms  need  to 
be  pre-processed,  including  DC  removal,  upsampling,  and  zero  padding,  to  minimize  frequency  leakage 
and  prevent  Fourier  transform  data  wrapping  during  the  dispersion  compensation  algorithm. 

Figure  5  shows  the  dispersion  compensation  and  dispersion  removal  results  of  a  350  kHz  SO  mode 
wave  on  a  3 -mm  aluminum  plate  after  propagation  distance  x  of  300  mm.  The  experimental  results  are 
found  to  be  very  close  to  our  numerical  prediction.  The  dispersed  SO  mode  wave  packets  were  well 
recovered  by  both  algorithms. 

Figure  6  shows  the  dispersion  compensation  and  removal  results  of  a  36  kHzv4o  mode  wave  on  a  1- 
mm  aluminum  plate  after  propagation  distance  x  of  400  mm.  The  experimental  results  are  very  close  to 
simulation  results.  Before  applying  dispersion  compensation  and  removal  algorithms,  the  AO  mode  wave 
packet  time  span  is  around  200  jus  .  In  contrast,  the  AO  mode  wave  packet  was  compressed  to  around  80 
jus  span  after  dispersion  compensation.  The  spatial  resolution  of  AO  mode  wave  packet  is  largely 
increased  after  applying  dispersion  compensation  and  removal  methods.  Again,  a  high  frequency 
component  was  observed  in  the  recovered  AO  when  using  dispersion  compensation  algorithm. 


^  ^  ^  1 

1  3. 5 -count  Tone  burst  at  350  kHz 

T 

V 

_ ^ ^ ^ _ 1 

0  30  60  90  120 


(a)  time,  us 


0  300  600  900 


time,  us 

Figure  5  Experimental  results  of  dispersion  compensation  of  350  kHz  SO  mode  on  a  3-mm  aluminum 
plate:  (a)  3.5-count  Hanning  windowed  tone  burst  center  at  350  kHz;  (b)  dispersed  SO  mode 
wave  after  x  =  300  mm  propagation  distance,  simulated  by  Eq.  (2);  (c)  recovered  SO  mode 
wave  in  spatial  domain  by  dispersion  compensation  algorithm;  (d)  recovered  SO  mode  wave 
in  time  domain  by  dispersion  compensation  algorithm;  (e)  recovered  SO  mode  wave  in  time 
domain  by  dispersion  removal  algorithm 
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Figure  6  Experimental  resuits  of  dispersion  compensation  of  36  kHz  AO  mode  on  a  1-mm  aluminum 
plate:  (a)  3-count  Hanning  windowed  tone  burst  center  at  36  kHz;  (b)  dispersed  AO  mode 
wave  after  x  =  400  mm  propagation  distance,  simuiated  by  Eq.  (2);  (c)  recovered  AO  mode 
wave  in  spatial  domain  by  dispersion  compensation  algorithm;  (d)  recovered  AO  mode  wave 
in  time  domain  by  dispersion  compensation  algorithm;  (e)  recovered  AO  mode  wave  in  time 
domain  by  dispersion  removal  algorithm 


4  APPLICATION  OF  DISPERSION  COMPENSATION  TO  PWAS  PHASED  ARRAY 


PWAS  phased  array  uses  a  group  of  sensors  located  at  distinct  spatial  locations  and  focuses  the 
guided  wave  front  in  a  specific  direction  by  sequentially  firing  the  individual  elements  in  the  array  at 
different  time.  The  beam  steering  and  scanning  is  implemented  using  the  embedded  ultrasonic  structural 
radar  (EUSR)  methodology  in  combination  with  the  PWAS  guided  wave  technology  (Giurgiutiu  and  Bao, 
2002)  (Yu  and  Giurgiutiu,  2007).  One  requirement  for  the  success  of  PWAS  phased  array,  a  single  mode 
non-dispersive  Lamb  wave  needs  to  be  employed  (Santoni  et  al.  2007).  This  largely  limits  the  operation 
frequency  selection  of  the  phased  array  though  the  single  mode  situation  has  been  solved  by  using  the 
frequency  tuning  approach.  In  this  section,  we  explored  the  application  of  dispersion  compensation 
algorithm  to  1-D  linear  PWAS  array  EUSR  interrogation  using  tuned  dispersive  Lamb  waves. 

Recall  that  EUSR  works  in  pulse-echo  transducer  mode  and  round-robin  data  collecting  pattern. 
For  an  M-PWAS  array,  Af  sets  of  signal  data  need  to  be  collected.  Assume,  the  f"  PWAS  is  the 
transmitter  supplied  with  tone  burst  excitation  and  the PWAS  is  the  receiver.  The  received  wave  can  be 
predicted  by 


^  In 


do) 


/,7  =  0,1,...M-1 


(7) 


where  TB{co)  denotes  tone  burst  excitation  in  frequency  domain;  r  is  transducer  location  vector.  To 
simulate  a  M  =  8  PWAS  array,  64  groups  of  source  data  will  be  generated.  Each  is  saved  in  .csv  data 
format.  For  EUSR  without  dispersion  compensation,  these  64  .csv  files  will  be  read  directly  and  imaged 
by  EUSR.  In  contrast,  for  EUSR  with  dispersion  compensation,  the  64  .csv  files  will  first  be  processed  by 
dispersion  compensation  algorithm  and  then  imaged  by  EUSR.  Four  cases  of  1-D  PWAS  phase  array  for 
damage  detection  were  simulated,  as  shown  in  Table  4.1.  Figure  7  shows  the  case  I-IV  aluminum  plate 
and  crack  setup. 


Table  4.1  Simulated  cases  of  damage  detection  with  1-D  PWAS  phase  array 


Lamb  wave 
freq.  &  mode 

Group 

velocity 

(m/s) 

Structure 

(Plate) 

thickness 

Crack  location(s) 
y crack  (mnVdeg) 

Sensor 

spacing 

(mm) 

Case  I 

3 -count  36  kHz  AO 

1162 

l-mm  Al 

200/90° 

8 

Case  II 

3.5-count  350  kHz  SO 

5242 

3 -mm  Al 

150/90° 

7 

Case  III 

3.5-count  350  kHz  SO 

5242 

3 -mm  Al 

150/90°,  170/90° 

7 

Case  IV 

3.5-count  350  kHz  SO 

5242 

3 -mm  Al 

150/90°,  166/90° 

7 

Figure  7 


1-D  PWAS  array,  crack  and  plate  setup:  (a)  case  I,  II;  (b)  case  III,  IV 


Figure  8  shows  the  EUSR  scanning  images  and  90°  A-scan  signals  of  one  broadside  crack  in  case  I, 
located  at  200  mm/90°  and  detected  with  36  kHz  AO  mode.  Without  dispersion  compensation,  the  EUSR 
scanning  image  is  big  and  blur  (Figure  8a)  because  its  A-scan  signal  at  90°  spreads  out  in  time  domain 


(Figure  8c).  In  contrast,  after  dispersion  compensation,  the  scanning  image  becomes  sharper  (Figure  8b) 
and  its  A-scan  signal  at  90°  is  compressed  (Figure  8d).  Similar  improvements  in  EUSR  scanning  images 
were  observed  when  using  350  kHz  SO  mode  to  detect  a  crack  in  case  II  with  the  aid  of  dispersion 
compensation  (Figure  9). 

Figure  10  demonstrates  the  EUSR  inspection  results  of  two  broadside  cracks  (case  III)  distanced  by 
20  mm  using  350  kHz  SO  mode  with  and  without  dispersion  compensation.  In  both  cases,  these  two 
cracks  were  detected.  However  a  higher  contrast  ratio  of  scanning  image  (Figure  10b)  was  obtained  after 
using  dispersion  compensation.  When  the  two  cracks  on  plate  are  separated  only  by  16  mm  (case  IV), 
EUSR  lost  its  resolution  without  dispersion  compensation,  i.e.,  only  one  crack  was  observed  in  canning 
image  (Figure  11a).  However,  after  applying  dispersion  compensation,  both  of  the  cracks  were  detected 
(Figure  11b). 
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Figure  8  EUSR  inspection  results  for  case  I  (3-count  36kHz  AO  mode,  crack  located  at  x  =200mm)  : 

(a)  mapped  EUSR  scanning  image  without  dispersion  compensation;  (b)  mapped  EUSR 
scanning  image  with  dispersion  compensation;  (c)  selected  A-scan  at  90°  without  dispersion 
compensation;  (d)  selected  A-scan  at  90°  with  dispersion  compensation. 


Figure  9 


EUSR  inspection  results  for  case  II  (3.5-count  350  kHz  SO  mode,  crack  located  at  x 
=150mm) :  (a)  mapped  EUSR  scanning  image  without  dispersion  compensation;  (b)  mapped 
EUSR  scanning  image  with  dispersion  compensation;  (c)  selected  A-scan  at  90°  without 
dispersion  compensation;  (d)  selected  A-scan  at  90°  with  dispersion  compensation. 


Figure  10 


EUSR  inspection  results  for  case  III  (3.5-count  350  kHz  SO  mode,  two  cracks  located  at  x 
=150mm  and  x  =  170mm)  :  (a)  mapped  EUSR  scanning  image  without  dispersion 
compensation:  (b)  mapped  EUSR  scanning  image  with  dispersion  compensation;  (c) 
selected  A-scan  at  90°  without  dispersion  compensation;  (d)  selected  A-scan  at  90°  with 
dispersion  compensation. 


Figure  1 1  EUSR  inspection  results  for  case  IV  (3.5-count  350  kHz  SO  mode,  two  cracks  located  at  x 
=150mm  and  x  =  166mm);  (a)  mapped  EUSR  scanning  image  without  dispersion 
compensation;  (b)  mapped  EUSR  scanning  image  with  dispersion  compensation:  (c) 
selected  A-scan  at  90°  without  dispersion  compensation:  (d)  selected  A-scan  at  90°  with 
dispersion  compensation. 

5  CONCLUSIONS 

Nonlinear  characteristic  of  the  wavenumber  is  the  cause  of  the  Lamb  wave  dispersion.  The 
dispersion  causes  the  elongation  of  received  waves,  deteriorates  the  spatial  resolution  of  the  waves,  makes 
the  experimental  data  hard  to  be  interpreted,  and  limits  the  selection  of  Lamb  wave  operation  frequency. 

In  this  paper,  dispersion  compensation  and  removal  algorithms  were  first  theoretically  investigated 
and  compared  for  recovering  SO  and  AO  mode  Lamb  waves  using  both  simulation  and  experimental  data. 
It  was  found  that  both  algorithms  were  able  to  well  recover  the  original  shape  of  a  dispersed  SO  wave 
packet.  By  mapping  a  wave  packet  from  time  domain  to  spatial  domain,  the  dispersion  compensation 
algorithm  is  also  able  to  directly  recover  a  wave  packet’s  spatial  location,  which  is  not  available  in  the 
dispersion  removal  algorithm.  However,  the  dispersion  removal  algorithm  outperformed  the  dispersion 
compensation  algorithm  for  recovering  the  AO  mode  Lamb  wave.  Moreover,  the  dispersion  removal 


algorithm  takes  less  computation  efforts  than  the  dispersion  compensation  algorithm.  Next,  the  dispersion 
compensation  algorithm  was  applied  to  EUSR  PWAS  phased  array  to  image  a  single  crack  and  two 
closely  located  cracks  on  a  plate.  From  the  simulation,  we  conclude  that  the  dispersion  compensation  can 
help  improve  the  EUSR  spatial  resolution.  In  the  presented  EUSR  dispersion  compensation  simulation, 
the  compensation  algorithm  was  applied  to  transducers’  raw  data  directly,  application  of  the  algorithm  to 
the  A-scan  data  is  suggested  for  future  work  to  save  computation  effort. 
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